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We present a quantum Boltzmann equation analysis of the spin-Hall effect in a diffusive Rashba 
two-dimensional electron system. Within the framework of the self-consistent Born approximation, 
we consider the roles of disorder-induced quasiclassical relaxation, collisional broadening of the 
quasiparticles, and the intracollisional field effect in regard to spin-Hall dynamics. We present an 
analytical proof that the spin-Hall current vanishes, independently of the coupling strength, of the 
quasiparticle broadening, of temperature and of the specific form of the isotropic scattering potential. 
A sum relation of the collision terms in a helicity basis is also examined. 

PACS numbers: 72.10.-d, 72.25.Dc, 73.50.Bk 

I. INTRODUCTION 

The spin-Hall effect refers to the appearance of a net polarized spin flow along the direction perpendicular to an 
external applied dc electric field. This phenomenon suggests a possible mechanism to control spin dynamics using 
an electric field and may be employed to resolve the challenge of spin injection in the emerging field of spintronics. 
In early studies, the spin-Hall effect was associated with the disorder-induced spin-orbit (SO) interaction, of an 
extrinsic nature—^ More recently, another disorder-free spin-Hall effect was predicted, respectively, in bulk p-doped 
semiconductors^ and in Rashba two-dimensional (2D) systems— This intrinsic effect essentially arises from disorder- 
independent SO interactions, such as the Rashba SO coupling, the SO interaction involved in the Luttinger model, 
etc. Experimentally, the spin-Hall effect has been observed in a bulk n-type semiconductor— and in a 2D heavy-hole 
system— 

In Rashba 2D electron systems, the spin-Hall current can be strongly affected by a spin-conserving interaction 
between electrons and impurities. This issue has been investigated in the diffusive regime extensively by various 
methods, including the Kubo formula^aiiifiiiiii* the spin-density matrix method^ and a nonequilibrium Green's 
function approach— Physically, to examine the effect of disorder on spin-Hall current, we must simultaneously 
take into account collisional broadening and quasiclassical relaxation, as well as the intracollisional field effectiii^i^ 
produced by the action of the electric field during the course of electron-impurity scattering. The most important 
disorder effect on spin-Hall current is associated with quasiclassical relaxation. This scattering process yields an 
additional term to the spin-Hall current, causing it to vanish— iiMiiiSiiiiii^ In addition, the collision- induced spread 
of the quasiparticle density of states (DOS)-namely, collisional broadening-also influences the spin-Hall effect. It 
reduces the value of the intrinsic spin-Hall conductivity with increasing broadening of the DOS.i^ For short-range 
disorder, the total spin-Hall current still vanishes when collisional broadening in both the intrinsic and quasiclassical 
terms of the spin-Hall current is considered— iiMLiS There is yet another disorder-induced quantum effect-namely, 
the intracollisional field effect-which can not be neglected— Li2*i2i During the scattering process, the DOS of the 
quasiparticles can be modified by the external electric field, whereupon relaxation becomes dependent on this field 
even in linear response. To date, only Sugimoto et al. have tried to investigate the spin-Hall effect including these 
three processes simultaneously— £ Performing a numerical calculation of the Keldysh nonequilibrium Green's function 
in a spin basis, they found a generally nonvanishing spin-Hall current: the value of the spin-Hall conductivity they 
obtained depends on the spin-orbit coupling constant, the transport lifetime, and the Fermi momentum. 

In this paper, we analytically investigate the spin-Hall effect in Rashba 2D systems using the nonequilibrium Green's 



function method, essentially the same as that used by Sugimoto et al.f& but in a helicity basis. We carefully analyze 
the nonequilibrium retarded and lesser Green's functions by simultaneously taking account of quasiclassical relaxation, 
collisional broadening, and the intracollisional field effect. We prove analytically that the spin-Hall current vanishes, 
irrespective of the spin-orbit coupling strength, of the form of the isotropic impurity potential, of temperature, and 
of the collisional broadening of the DOS. We also discuss a sum relation of the relaxation terms in the helicity basis. 

This paper is organized as follows. In Sec. II, we present the kinetic equation for the "lesser" Green's function. 
The vanishing of spin-Hall current is proven in Sec. III. Finally, we review our results in Sec. IV. Several appendixes 
provide requisite details which could be read before Eq. J3J by readers desiring a "proof first" exposition. 
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II. KINETIC EQUATION FOR A 2D RASHBA SYSTEM 

We consider a 2D Rashba electron system with a single-particle noninteracting Hamiltonian given by^ 

P 2 

h a (p) = — + ap-(nxa). (1) 

Zm 

Here, p = (p x ,p y ) = (pcos^> p ,psin^> p ) is the 2D electron momentum, m is the effective mass, a = (o~ x , o~ v , a z ) are 
the Pauli matrices, a is a spin-orbit coupling constant, and n is a unit vector perpendicular to the 2D electron plane. 
By a local unitary spinor transformation U pi 

U = — ( 1 1 

the Hamiltonian Q can be diagonalized as fto(p) = £/p t "/io(p)^p = diag(ei(p), £a(p)) with e M (p) = |^ + (— l) M ap 
(/i = 1, 2) as dispersion relations of two spin-orbit-coupled bands. In this paper, we investigate the spin-Hall effect in 
a 2D Rashba electron system driven by a constant, uniform electric field E along the x axis. In the Coulomb gauge, 
the electric field can be described by a scalar potential, V = — eE ■ r, with r as the electron coordinate. 

In a realistic 2D system, electrons experience scattering by impurities. We assume that the interaction between 
electrons and impurities can be characterized by an isotropic potential, V(|p — k|), corresponding to scattering of an 
electron from state p to state k. 

The nonequilibrium Green's functions G p !< in spin basis for 2D electron systems with Rashba SO coupling are 
defined as usual. They are 2x2 matrices and obey Dyson equations presented in Appendix A. For brevity, hereafter, 
we employ a subscript p to denote the arguments of the Green's functions and self-energies, (p, to). It is most 
convenient to study the lesser Green's function in the helicity basis, Gp = U p G p U p , in which the unperturbed lesser 

Green's function is diagonal. Using the transformation from spin basis to helicity basis, the Dyson equation for G p 
can be written as 




ieE- (J- +p ^ G< +ieB.[G<,V p U+U p }+ap[a z ,G<} + l -^eE- |-B< = (3) 

with Ap = [N p ,Gp], Bp = {N p ,Gp}, N p = U p (n x a)U p , and as a scattering term. In general, I\ has a 
complicated form because the momentum and time variables of the Green's functions and self-energies are intertwined 
due to the presence of the electric fields Physically, this feature involves the intracollisional field effec1*i£*i2*i2i and 
results in an additional electric-field-dependent scattering. In this paper, we restrict our considerations to the linear 
response regime. Based on this, the scattering term can be expressed as the sum of three terms: If = In + + lis- 
The first term, In does not involve the explicit appearance of the electric field (although it is implicit in G^ and self- 
energies S^). The terms J/2 and I13 arise from the first-order gradient expansion involving explicit linear dependence 
on the electric field, and they are a manifestation of the intracollisional field effect. The explicit forms of quantities 
In and I12 are similar to In and I12, respectively, but with helicity-basis Green's functions and self-energies in place 
of the corresponding spin-basis ones (the forms of In and In are provided in Appendix A). The term I13 arises from 
the local nature of the transformation and its form is obtained from I12 by using | V P P [- • • , o~ x ] instead of operator 

v,». " ' " 

From the explicit form of the scattering term, it may be seen that Ii relates to the retarded Green's function G p in 
which the collisional broadening is considered. Detailed properties of G p are analyzed in Appendix B. We find that 
the linear electric field correction to the retarded Green's function G\ p in a 2D Rashba system does not vanish in 
general (although it has vanishing diagonal elements). Also, we make clear that disorder scattering is not essential 
for the nonvanishing of G^ p : G^ p exists even when the electron- impurity collisions are ignored. 

It is well known that the linear electric field correction to the retarded Green's function vanishes for a one-band 
electron gas^i However, for two-band systems, if a transition between the two bands is present, G\ p no longer vanishes 
(see, for example, Ref.l22l p. 215). In the 2D Rashba system that we study, there are two spin-orbit-coupled bands 
arising from the structure asymmetry and a polarization process can be induced between these two bands by a dc 
fields Accordingly, the linear dc electric field correction to G p in a 2D Rashba electron system is nonvanishing. 

In Eq. 10, the electron-impurity scattering is embedded in the self-energies Ep^. In the present paper, we consider 
the disorder collisions only in the self-consistent Born approximation. It is widely accepted that this treatment is 
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sufficiently accurate to analyze transport properties in the diffusive regime. On this basis, the self-energies in helicity 
basis take the forms 

E^<=n i ^f(p,k)G^ < f+(p ) k), (4) 

k 

where T(p,k) = U + (k)V(p — k)f/(p) and rij is the impurity density. 

Further, we assume that the applied dc field is weak enough that only the linear response is required to be considered. 
After linearizing Eq. , we can write the kinetic equation for the linear electric field correction to the lesser Green's 
function, Gf p = G p — Gq p , as 

- apCx + ieE • V p G< p - ^eE ■ V p ^ p A, + R = I^\ (5) 



where the matrices C± and Dq are given by 



D Q = f ( G 0p)ll _ ( G 0p)22 \ ,y> 

\ ( G 0p)22 - (Gop)ll J 

R is a remainder term which can be expressed as a total derivative with respect to ui. Gq p is the unperturbed 
equilibrium lesser Green's function. It is a diagonal matrix and simply relates to the retarded Gq p by the Kubo- 
Martin-Schwinger relation 2 ^ 

G< = -2m F (a;)ImGS p , (8) 

with tif(w) as the Fermi function. 

The scattering term ij 1 ' is the sum of , 1^ , and 1^' , the linear electric field parts of In, /; 2 and I13. The 

explicit forms of these quantities are presented in Appendix C. It should be noted that J.g and Jjg are diagonal 
and off-diagonal matrices, respectively. In deriving Eq. J3J), the properties of the retarded Green's function in the 
presence of collisional broadening are employed: its linear electric field correction is an off-diagonal matrix, while the 
unperturbed one is diagonal. 

From Eq. (J5J we can formally express the off-diagonal element of G^, as 

(Gf p )i 2 - ^{^E.V p p [(G o < p )h-(G o < p ) 22 ] 

+(/, (1) ) 12 -i? 12 }. (9) 

At the same time, the diagonal elements of D 1 ' take the form, 

= icE ■ V p (G< p ) w + R^, (10) 

with /i = 1, 2. 

III. VANISHING SPIN-HALL CURRENT 

We are interested in a spin current polarized along the z-direction and flowing along the y axis when a dc electric 
field E is applied along the x axis-ie. J*. In spin basis, it relates to the single-particle spin current operator j* by 
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There has been some debate about the definition of the spin current j* in recent years £££§£±2^ It seems reasonable to 
define the spin current from a continuity equation relating to the spin, as shown in Ref.l28t Conventionally, the single- 
particle spin current in spin basis is defined as j* = \{v yi a z } with the velocity operator v y . Zhang et al. suggested that 
the spin current operator has an additional term arising from the spin torque density, (jy) a dd = (j//2i)[cr z , /io(p)]- 28 
Actually, the contribution from this term to spin-Hall current vanishes in a homogeneous 2D Rashba electron system 
to linear order in the dc field. This can be seen from the fact that the additional spin-Hall current operator, explicitly 
given by (jy) a dd = ctyP 1 ® — ~^ a "^~P ' a : 1S invariant under inversion of momentum p, i.e. p — > —p. However, p 

reversal leads to a sign change of the linear electric field correction to the lesser Green's function: Gf_ p = — Gf 
[this result can be seen from Eq. JSJ by inverting the momentum]. Based on this, it follows from the p integration of 
Eq. Ijlljl that the contribution to spin-Hall current arising from the additional spin-current term vanishes. Hence, we 
only need to consider the conventional component of the spin-Hall current operator. 

Consequently, we write the spin current operator J* in the simple helicity-basis description as 15 

^=E/^|MG P <),, (12) 

P J 

Substituting Eq. J^J into this equation, we obtain 



J y=E 



du> 1 

— sin0 p eE • V p 0plm 

Z7r 4ma 



(G 0p )ii - (G 0p )22 | + 



E 



du 1 
27r 2ma 



sin </>pIm(/ i (:i '') 1 2 



(13) 



In this, the u> integration over the term R vanishes. 

To further simplify the expression for J*, we analyze the second term on the right-hand side of Eq. In 

connection with this, we find that there is a sum relation among the elements of 2, in the self-consistent Born 
approximation: 



£/isin^Im(/P) 12 =X; 



du cos 4 



Re(4 1) ) 11 -Re(4 1, )22 



(14) 



Furthermore, separate relations also hold for each component of J ; 



(i). 



2tt 



clio cos 4 
2 



Re(/( L 1) ) 11 -Re(/^) 2 2 



(ih 



(15) 



J2 1 ^z sin( t ) p lm ( I i3 ) )^ = J2 



2tt 



du: cos < 



Re(4 1) ) 11 -Re(4 1) ) 2 2 



The detailed proof of these equations is presented in Appendix D. 



Based on the sum relation Eqs. I|14fl and (|10fl for diagonal components of D X ' , we obtain 

dm sin < 



J y = E 



2. 4ma- eE - V ^ Im 



E 
p 

E 

p 

E 



dw cos i 



27r 4ma 



(Gop)ll - (Gq p )22 

^eE • (Vplm [(G< p )h - (G< p ) 22 



dio eE • V p cos <A P 
-lm 

dwcos^ p / 

-eE • Vnlm 



2ir Ama 



(Gop)ll - (Gq p )22 



(Gop)ll - (Gq p )22 



(16) 



(17) 



It is apparent from Eq. 117|) that the integrand on its right-hand side becomes a total derivative with respect to p. 
Consequently, the vanishing of the spin-Hall current is obtained immediately upon the momentum integration. We 
note that this analytical result is significantly different from that of the numerical analysis of Sugimato et al.fi& which 
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purports to exhibit a nonvanishing spin-Hall current. Our analysis shows that this is not the case even with the "new" 
spin current term included. 

It is obvious from Eq. (|17|l that the spin-Hall current vanishes independently of the strength of SO coupling, of 
the specific form of the isotropic scattering potential, and of the lattice temperature. Moreover, it is valid even 
for arbitrary broadening of the DOS. However, this result holds only in the diffusive regimes-namely, in the regime 
lokp > 1 (Id is diffusion length and kp is the Fermi wave vector), because we restricted treatment of the electron- 
impurity scattering to the self-consistent Born approximation (otherwise, contributions from the maximally crossed 
diagrams of electron- impurity scattering would have to be taken into account). 



The spin-Hall effect in 2D diffusive Rashba systems has been investigated by means of a quantum Boltzmann 
approach in this paper. In the self-consistent Born approximation, we have considered all effects induced by electron- 
impurity scattering-namely, quasiclassical relaxation, collisional broadening, and the intracollisional field effect. We 
have proven analytically that the spin-Hall current vanishes, irrespective of the strength of SO coupling, of the 
broadening of DOS, and of the form of the isotropic scattering potential. The sum relation of the scattering terms in 
the helicity basis has been analyzed. 

Our treatment is applicable at nonvanishing temperatures (not restricted to T = 0), and it is not restricted in 
anyway to short-range scattering potentials. We have also made it clear that the spin-Hall effect vanishes even when 
correctly considering the influence of the electric field in the scattering process (intracollisional field effect). 

It should be noted that this proof of vanishing of spin-Hall current also applies for a Dresselhaus linear-momentum 
SO coupling, since it relates simply to the Rashba SO interaction by a global transformation. Furthermore, we have 
also examined the spin-Hall effect when both the Rashba and Dresselhaus SO interactions are active, again finding 
zero spin-Hall current. 

Nevertheless, we note that our proof of vanishing of the spin-Hall effect has its validity limited to 2D systems with 
linear-momentum SO coupling. Also, this proof requires ho to be parabolic in the absence of SO coupling. Otherwise, 
the spin-Hall effect should be examined anew. This is to say that our result does not conflict with the nonzero 
experimental spin-Hall observations^: the Hamiltonians of bulk and 2D heavy-hole systems in which nonvanishing 
spin-Hall currents were observed differ significantly from the Hamiltonian ft-o(p) studied in this paper. 

Note added. After this work was completed and submitted, an unpublished paper of Krotkov and Das Sarma 
reported a nonvanishing spin-Hall effect in 2D Rashba electron systems by considering a nonparabolic effect^ 
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APPENDIX A: NONEQUILIBRIUM GREEN'S FUNCTIONS IN SPIN BASIS 



In spin space, the 2x2 nonequilibrium Green's functions G(ri, n; T2, ts) obey Dyson equations given by2^ 





(Al) 
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6(1 - 2) + J dr' J dr'G(n, n; r', r')£(r', r'; r 2 , r 2 ). 



(A2) 



In these equations, the electron-impurity interaction is embedded in the self-energies, £(r', r'; r 2 , T2). 

To further simplify Eqs. (|A1|) and (|A2|) . we introduce relative and center-of-mass variables r = ri — r 2 , t = t\ — r 2 , 
R = (ri + i*2)/2, and T = (ri + Ta)/2. Following this, we construct gauge-invariant retarded and lesser Green's 
functions G^ in momentum-frequency space following established procedures^ Under homogeneous and steady- 
state conditions, these Green's functions satisfy the following equations: 



P 

2m 



G p 



ip B' 



Irl 



«r2, 



(A3) 



eE ■ 



duj 



ip - A < 



In 



(A4) 



where A p < 



[n x a, G^] and B p !< = {n x <r, G^} arise from the Rashba spin-orbit interaction term in the 
Hamiltonian. The first terms I r \ and In on the right-hand sides of Eqs. I|A3(1 and i|A4|l . respectively, do not involve 
the explicit appearance of the electric field 



Irl — ^p^p 



G~ r v» r 
P P' 



(A5) 



in — ^p^p — b p z, p — bp^p -t- z^p bp. 



(A6) 



The terms I r 2 and 7/ 2 arise from the first-order gradient expansion involving explicit linear dependence on the electric 
field, 



d A , 



d 



Jra = -eE • V n E^ — G^ - — S^V n G: 



p^p 



<9w 



<9w 



j p v p^p 



(A7) 



i/2 



-eE- 



|6;v P E< + v P E<|-e« 



(A8) 



Incidentally, ignoring the collisional broadening, i.e., the terms involving d/duj, and employing the generalized 
Kadanoff-Baym ansatz, Eq. i|A3|l reduces to the kinetic equations of the distribution functions presented in the 
previous studies^ 



APPENDIX B: RETARDED GREEN'S FUNCTION 

To analyze the properties of the retarded Green's function, we begin from the Dyson equation in the hclicity basis: 



V 

2m 



G p + apK,G p } + y eE-— A p 



2 + Irl + Ir 



(Bl) 



This equation is obtained from Eq. I|A3(1 by the transformation from spin basis to the helicity basis. The quantities 
I r i and I r 2 are the helicity-basis analogues of the quantities I r i and I r 2, respectively. I r 3 has a form similar to /,. 2 , 
but with §V p 0p[- • • , a x ) in place of the operator V p . 
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Obviously, in the helicity basis, the noninteracting retarded Green's function is diagonal 

9p = diag ((uj - Ei(p) + iS)- 1 , (u - g 2 (p) + iSy 1 ) . 



(B2) 



Including collisional broadening in the definition of the unperturbed retarded Green's function Gg p , it obeys the 
Dyson equation in the absence of the electric field, as 



V 

2m 



g; 



ap 



1 



up t I^Wzi Gq p } — 1 + ^(^opGSp + GgpSop) 



(B3) 



The solution of this equation is also a diagonal matrix, (Gq p )i2 = (Gqp)2i ~ 0> with elements independent of the 
direction of momentum. To verify this, we examine the form of the self-energy Sq p below. 

Let us first discuss the forms of self-energies tZ; < in general. Installing the actual form of matrix U p into Eq. 10}, 



we find 



r,< 
k 



W3[0'iiGi < ] 



(B4) 



Explicitly, the self-energies can also be rewritten as 



(^p )/w 



^inp-k)| 2 { ai (G^U+ 

k 

a2(G£ < ) jiZj ? + ^[(G^)^ - (G^)^]} 



(B5) 



In these expressions, fi = 3 — ji and <Zj(« = 1,2,3) are factors associated with the directions of momenta, a\ = 
1 + cos(<j!>p - fa), a 2 = 1 - cos(0 p - k ), and a 3 = sin(0 p - k ). 

From Eq. (|B5|) . we see that the diagonal elements of the Green's functions appearing in nondiagonal elements of 
S p !< are always combined with the factor 0,3 . If the solution of Eq. (|B3|) . GJj p , is diagonal and independent of 
the direction of momentum, the self-energy term involving 03 vanishes under momentum integration. Hence, the 
self-energy, Sg p , also becomes diagonal and independent of momentum direction. As a result, Eq. (|B3|) reduces to a 

momentum-direction-independent equation for the diagonal elements of Gg p . This implies that the assumed form of 

Gq p is consistent with Eq. i|B3|) . From uniqueness of the solution of the Dyson equation we immediately conclude that 

such a Gq p is just the needed solution. It should be noted that, physically, the momentum-direction independence of 
the unperturbed retarded Green's function is associated with the rotational symmetry of the unperturbed Hamiltonian 
in helicity basis. 

Considering the angular independence of Gq p , an important relation can be derived: 



£ \V(P - k)| 2 sin(0 p - k )[(GS k ) 22 - (Ggjju] = 0. 



(B6) 



This equation is needed to prove the sum relation of the scattering term (Appendix D). 

To the first order in the dc field, the Dyson equation for the linear electric field correction to the retarded Green's 
function, G£ p = G p — Gg p , can be written as 



2m 



(gDm, + (i - <W 



d r 



(Gq)h - (Go)22 



aeE 



Ify, Ir2 t an d are * ne linear electric field parts of J r i, I r 2 and I r ^. 1^ has the form 



f(i) 

'rl 



(Go p )ii(Si p )ii^ 
[(Ggp)n + (G5 P ) 2 2](£i p )2i 



[(Gop)n + (G5p)22](SJ p )i2 

(0^)22(^^)22 



(G 



Op 



(B7) 



(B8) 
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where SJ p = ^ k T(p, k)G^ k T + (p, k). 1$ vanishes due to the diagonal nature of the unperturbed retarded 
Green's function and self-energy, f^g is off-diagonal and symmetric, (1^)12 = (-^3^)21- From this Dyson equation 
and the explicit forms of self-energies l|B5(l . we find that G\ p becomes an off-diagonal matrix with symmetric elements, 
(G' lp )i2 = (Gj p )2i- Also, E£ and hence I^i have vanishing diagonal elements. At the same time, Eq. (|B7|I reduces 
to an equation for nondiagonal elements of G^ p . 

Actually, G^ p depends on the direction of momentum through a sine function. The expression for G^ p can be 
formally written as 



(^ip)w — 



eE sin < 



— p 2 /2m — t p 







~2 oTjK ^o P )ii ~ (Gq p )22] + [(Go p )ii + (Gg p ) 22 ]A p 



(B9) 



with A p = 2n^ k |y(p - k)| 2 cos(0 p - k )(^ p )i2, r" 



= 2n t j: k \V(p - k)| 2 [(G5 p ) n + (G5 p ) 22 ], and Q\ p = 

Gi p /[eE sin(f)p] as a quantity independent of the direction of momentum. In the case of short-range disorders, 

the A p and t p become momentum-independent constants. From Eq. (|B9|I it is obvious that G\ p does not vanish. 

We note that the electron-impurity scattering is not essential for the nonvanishing of G\" p . In this, the second 

term on the left-hand side of Eq. (|B7|I plays a key role. In the absence of the electron-impurity collisions, G\ p can be 
obtained analytically as 



(G' lp )i2 — (G[ p ) 2 i — - 



aeE 



A{u — p 2 /2m] 



■ sin 



h 9 



(Go P )n - (Gq p ) 2 2 



(BIO) 



Incidentally, from Eq. (|B10|) we can derive the retarded Green's function presented in our previous work^ Using the 
equality (Gq p ) = (gL,) in the absence of electron-impurity collisions, we rewrite Eq. (|B10|> as (GJ p )i 2 = (G^ p ) 2 i = 
Ai + A 2 with 



and 



A 2 



eE 
Aap 2 



eE . L d 



(G 



Op) 



(G 



Op) 



(Gq p )h + (G 



Op) 



(Bll) 



(B12) 



Ai is just the term that we employed in Ref.ll5l and A 2 is a total derivative with respect to lo. In the absence of 
electron-impurity scattering, the term A 2 has no effect on the spin dynamics of the considered system, because the 
quantity directly associated with spin-Hall current, Gf p , takes the form 



-2inF(u>)IuiAi + TZ, 



(B13) 



with TZ as a remainder term which can be expressed as a total derivative with respect to to [this result is obtained 
from Eq. JjJJl]. 



APPENDIX C: EXPLICIT EXPRESSION FOR 

The scattering term 1^ can be expressed as the sums of ij^ , In i ano - ^13^ the corresponding linear electric 
field parts of In, 1^, and 7; 3 . Further, ijP can be rewritten as iffl — F% + F 2 , with an off-diagonal matrix Fi = 
K^o P - G<S? p - 6j p £< + E<G? p , or explicitly, 

(Fi)fj.p. — (Sip)/ip(Go p )pp — (E" p ) A( p(G^ p ) AtAt — (Gi p ) ^(TiQp) + (G" p ) A1 p(EQ p ) A1A1 , (CI) 

and a general matrix F 2 with elements 

= 2i[Im(ES p ) w (G< p ) w - Im(GS p ) w (E< p )p M ], (C2) 
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where linear electric field correction to the lesser self-energy, E< p , takes the form: E< p = ^ k T(p, k)G^ k T + (p, k). 
Furthermore, the diagonal jig and off-diagonal Z,g matrices have the forms 

(IlP)w = «eE • f VpRe(So p ) w — (G^ p ) w - — Rc(So p ) w V p (G5p) MM 

-V p Re(GS p ) w £(E <) w + A Re( G 5p ) w Vp(E< p ) w ) , (C4) 



(Gop) w ^jR< e (Sop)i^ — ^j(Go P )^R e (^opW 



^jRc(G 0p ),yy (EQp)^^ + Re(G 0p ) I yy^j(I]Q p ) ! y !/ 



(&W> (C5) 



From the definition of the lesser Green's function, we know that (Gf p ) M „ = —(Gf p )l^- This implies that its 
diagonal elements are pure imaginary- i.e., Re(Gf p ) MM = 0. Hence, we see from Eq. (|B5|) that the diagonal elements 
of the linear electric field part of the self-energy, E^, are also pure imaginary. Consequently, the diagonal matrix F 2 
is pure imaginary, as well as the matrix iffl- 

APPENDIX D: PROOF OF SUM RELATION 

First we prove Eq. (|15(l . Since Im(f/ 1 2 ' ) )i 2 is the sum of (Fi)i2 and Im(i*2)i2, we first consider the integral X\ = 
2p sm< Ap(-^i)i2- Explicitly, it takes the form 

l x = ]Tsin0 p {Re(^ p ) 12 [Im(G< p ) 22 - Im(G<) n ] - Re(GJ p ) 12 [Im(S < p ) 22 - Im(E < ) n ]} . (Dl) 
p 

Substituting the self-energies into Eq. (|D1|I . we have 

T x = ^sin0 p |y( P - k)| 2 {Re(G^ k ) 12 [Im(G< p ) 22 - Im(6$,)ii] 

6 P - ^ k )Re(G^ p ) 12 [Im(G< k ) 11 - Im(G <) 22 ]} . (D2) 

Interchanging the dummy integration variables p and k, X\ can be rewritten as 

Ii = \ V (\P - k l)| 2 ( sin ^P ~ sin0 k cos(0 p - k ))Re(G^ k ) 12 [Im(G< p ) 22 - Im(G <) u ] 
P k 

= l^dP - k ')| 2 cos0 k sin(0 p - k )Re(G^ k ) 12 [Im(G< p ) 22 - Im(G< p )n]. (D3) 

Pk 



pk 

cosf 



From relation i|B6() it can be seen that T\ — 0. 

Next, we analyze the remaining component of the integral, X 2 = J^ p sin0 p lm(i 7 2 ) 12 , 

1 2 = ^sm^ p {Re(G< p ) 12 [Im(E5 p ) 11+ Im(SS p ) 22 ] 
p 

+Im(G< p ) 12 [Re(ES p ) 11 + Re(E£ p ) 22 ] - (E <-> G)} . (D4) 
Again substituting the self-energies into Eq. l|L)4|i . this integral can be simplified as 

I2 = ^|F(p-k)| 2 cos0 k sin(0 p -0 k )Re(G< p )i 2 [Im(GS k )ii+Im(GS k ) 22 ] (D5) 

P,k 

-- ^[cos0 k - cos (f)p cos(</)p - </> k )][Im(Go p )ii + Im(Go p ) 22 ] [Im(Gf k )n - Im(Gf k ) 22 ]. 

P,k 
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Here, we have used Eq. I|B6(1 . In this, the vanishing of the real parts of the diagonal elements of the distribution 
function- i.e., Re(G^~ k ) w — 0-has also been considered. Combining the terms proportional to cos0 p cos(0 p — 0w_) in 
the second line of Eq. (|D5|I with a similar term in the first line, we obtain 



J 2 = ^cos0 p {[Im(S< p ) n -Im(E< p ) 22 ][Im(GS p )n+Im(GS p ) 22 ] 
p 

-[Im(S<)n +Im(E < p ) 22 ][Im(G^ p )n-I m (GIp) 22 ]} . 
Further, we note that in the self-consistent approximation, there is a vanishing quantity K. given by 
K = X>s^ p {[MGS P )n - Im(G5 p ) 22 ][Im(E< p ) n + Im(E< p ) 22 ] — (£<-> G)} . 



(D6) 



(D7) 



The vanishing of K, = can be shown by inserting the explicit forms of the self-energies, Eq. 1B5I) . into JC and 
employing Eq. (|B6|) . Finally, adding JC to the right hand side of Eq. (|D6|) , we find 



(D8) 



12 = ^ cos0 p (-iy +1 {lm(Sf p ) w Im(G5 p ) w - Im(£< p ) w Im(Gy MM ) 



cos 



£^[Re(F 2 ) u -Re(i* 



2)22\, 



which is Eq. (|T5|) . 

In the following, we prove Eq. I|ltj|) . From Eq. (|C5f) . the involved integral, denoted as I3, can be written as 



T 3 = / ^(-ircos0 p eE-V p {(G<) w ^-Re(ESW- 



aw (G <)„ 



Re(ES) 



— RetGoW 



(E^)^^ + Re(GQ) l / l /^j(S^ : ) !v;/ 



(D9) 



Here, the momentum integral of Eq. (|C5|I has been performed by parts integration, and the contribution from (i? 2 )i 2 
vanishes due to its total derivative nature with respect to lu. Each part of the integrand involves two functions Gq^ 
and Sg^. Further integration by parts with respect to to in terms of the form (^;V p A) B ([with A and B being 
(Gg^)^ and (Eq < ) vu , respectively, or vice versa] results in the operators V p and d/dw acting individually on A or 
B-i.e., (£jA) V P B. Thus, we obtain 



E 



did 

2^ 



(-l)"cos (j) p eE ■ 



V p (G<) w A Re (S5) w - A( G <) w V p Re(S5) 



-^Re(Gl) vv V v {t<U + V p Re(GS) w ^ (E<) w 



(D10) 



Using Eq. Eq. 

can be derived immediately. 
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